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^ ■ Abstract 

in 

C^^ \ In this note, we construct simple stochastic toy models for holographic gauge theo- 

1.^^ ' ries in which distributions of energy on a collection of sites evolve by a master equation 

■' ■ with some specified transition rates. We build in only energy conservation, locality, and 

Q ! the standai^d thermodynamic requirement that all states with a given energy are equally 

cn ' likely in equilibrium. In these models, we investigate the qualitative behavior of the 

dynamics of the energy distributions for different choices of the density of states for 
the individual sites. For typical field theory densities of states (log(p(£')) ~ £'°*^^), 
S^ \ the model gives diffusive behavior in which initially localized distributions of energy 

spread out relatively quickly. For large N gauge theories with gravitational duals, the 
density of states for a finite volume of field theory degrees of freedom typically includes 
a Hagedorn regime {\og{p{E)) ~ E). We find that this gives rise to a trapping of energy 
in subsets of degrees of freedom for pai^ametrically long time scales before the energy 
leaks away. We speculate that this Hagedorn trapping may be part of a holographic 
explanation for long-lived gravitational bound states (black holes) in gravitational theo- 
ries. 



^ 



1 Introduction 

According to the AdS/CFT correspondence in string theory [ID 13, gravitational dynamics 
can arise as an emergent phenomenon in the physics of strongly-coupled large A^ gauge 
theories. While there is now a great deal of evidence for the validity of the correspondence, 
the reason for the emergence of a dual spacetime and its associated gravitational physics is 
still not well understood. 

An intriguing aspect of the correspondence is that many different field theories in dif- 
ferent numbers of dimensions are believed to have gravity duals. The physics of these field 
theories can be quite different, and details of the dual gravitational theories vary, but a very 
large class of them include Einstein gravity as part of the dual spacetime physics. Thus, 
Einstein gravity (including the physics of black holes) seems to emerge as a universal sector 
in the physics of a large class of strongly-coupled large N field theories. 

A possible explanation for this universal emergence of gravity is that it is not directly 
associated with any detailed dynamical properties of the field theory, but rather to "ther- 
modynamic" effects, associated with the very large number of degrees of freedom present 
for large A^ theories and the large energies (of order A^^) required for any non-trivial bulk 
gravitational dynamicslll In this picture, the dual spacetime metric is thought of as a set of 
macroscopic thermodynamic functions associated with the field theory (similar to energy 
density or pressure), and the large number of degrees of freedom guarantees that the future 
evolution of these macroscopic variables (for typical states) is completely determined by the 
present values through a set of simple macroscopic equations, without further knowledge 
about the detailed microscopic state. 

Furthermore, the microscopic physics enters into the macroscopic equations only in a 
limited way; for example, in the heat equation, the microscopic physics determines the ther- 
mal conductivity as a function of temperature. Once this single function is determined from 
the underlying microphysics, the dynamics of heat flow in the system is determined. Sim- 
ilarly, it may be that relatively few details of the microscopic physics of large A^ gauge 
theories are enough to determine the macroscopic equations that govern the evolution of the 
metric in the dual spacetime. 

Motivated by these observations, we construct in this paper some very crude toy models 
for large A^ gauge theories, and ask whether we can reproduce some of the physics associated 
with the existence of a dual gravitational description by building in only a few key features 
of the original gauge theory. In particular, we consider stochastic models in which energy is 
distributed between degrees of freedom living at various sites, and allow these configurations 
to evolve randomly, demanding only that energy transfers are restricted to neighboring pairs 
of sites. Our model is classical and stochastic, but we emphasize that it can still incorpo- 
rate quantum effects through the choice of the density of states and the transition rates. By 
analogy with ordinary thermodynamics and hydrodynamics, we expect that by focusing on 
just the dynamics of the energy, a classical stochastic picture should suffice to describe the 



'Many connections between gravity and thermodynamics have been made in the past, both within the 
AdS/CFT correspondence and independent of it. In particular, the physics of black holes is governed by laws 
that are in perfect analogy with the laws of thermodynamics. A very incomplete list of references is ||3]@]|5]|6] 
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physics. In essence we are integrating out all degrees of freedom besides the energy den- 
sity and studying the reduced stochastic dynamics; quantum effects enter indirectly through 
plausible choices for transition rates and the density of states. 

We find that building in only a single feature related to large A^ gauge theories - a suf- 
ficiently rapid growth in the density of states for the degrees of freedom at a site - leads to 
a trapping of energy in localized regions for parametrically long times. These results sug- 
gest a direct connection between the Hagedom (or super-Hagedom) density of states found 
for certain subsets of degrees of freedom in large N gauge theories and the long lifetime of 
gravitational bound states (black holes)B The Hagedom density of states is also famous as 
the density of states associated with perturbative strings. Thus, we expect that any theory for 
which weakly coupled strings play a dominant role in some regime must have a Hagedom 
density of states for some subset of its degrees of freedom in some energy regime. Our toy 
models suggest that these theories will be able to trap energy for long time scales in subsets 
of the degrees of freedom, so from this perspective, it is natural that string theories are able 
to describe black holes. 

More generally, we find that for slower-than-Hagedom growth, energy diffuses; small 
fluctuations about a uniform energy distribution are govemed by an ordinary heat equation. 
For faster-than-Hagedom growth (found for a certain range of energy in theories whose 
dual gravitational theory includes small black holes), a clustering behavior is observed; for 
example, small fluctuations about a uniform energy distribution are govemed by the time- 
reversed heat equation. This clustering behavior for super-Hagedom growth in the density 
of states may be part of a field theory explanation for the formation of small black holes in 
the dual theory. 

In gauge theories whose gravity dual admits long-lived evaporating black hole states 
(e.g. in confining gauge theories where a bulk infrared wall prevents small black holes from 
falling through a Poincare horizon), the density of states for a localized subset of field theory 
degrees of freedom typically includes several regimes, with slower than Hagedom growth 
at the smallest and largest energies and an intermediate regime characterized by Hagedom 
or faster-than-Hagedom growth. Inserting such a density of states into our stochastic toy 
model, we find that a large amount of energy localized to a small region will initially spread 
out to some region of a characteristic size determined by the energy and then remain there 
for a long time, very slowly losing energy to the surrounding degrees of freedom. Thus, we 
find behavior qualitatively similar to the formation and evaporation of a black hole. 

A physical picture of our results is that the energy is getting trapped within some set of 
complicated local degrees of freedom that can only be accessed at sufficiently high energy. 
At these energies local intemal rearrangements of the degrees of freedom are much more 
likely than large scale spatial rearrangements and diffusion is inhibited. While ordinary 
conformal field theories and confining gauge theories may not have such a set of many local 
degrees of freedom, such a picture does emerge from large A^ confining gauge theories as 
we discuss in more detail below. Roughly speaking, once the energy density is above a value 
set by the confinement scale the system may access a large A^^ number of local degrees of 



^For a Hagedorn density of states, the density of states grows exponentially with energy as p{E) oc 
exp(%£;). 



freedom and the phenomenon of getting trapped in Hilbert space can occur. 

The connection between black holes and Hagedom or super-Hagedom densities of states 
is certainly not new. It is well known that the number of black hole states grows very rapidly 
with energy, and it has been suggested that this rapid growth may be directly related to the 
Hagedom growth of perturbative string states (9]. In the AdS/CFT context, the appearance 
of Hagedom densities of states in field theory has been linked to the existence of a deconfine- 
ment phase transition as temperature is increased; this deconfined phase of the field theory 
corresponds to large black holes or black branes in the dual gravitational theory (see, for ex- 
ample ifTOlfTTlfTll '). The novel aspect of the present work is that we relate the behavior of the 
density of states to the dynamics of energy transfer between degrees of freedom in the field 
theory. We are thus able to suggest a mechanism by which energy can be trapped for long 
periods of time in localized regions of an interacting, translationally-invariant open system, 
providing a possible field theory explanation for the existence of long-lived (but unstable) 
black hole states in the dual theory. 

The paper is organized as follows. In section 2, we motivate and describe the basic 
stochastic models that we consider. We show that Hagedom densities of states are naturally 
associated with static or nearly static energy distributions. In section 3, we derive differential 
equations that govem the evolution of the energy distribution in the thermodynamic and 
continuum limits of our model, and investigate how the growth rate in the density of states 
affects the evolution. In section 4, we review relevant results on the density of states for field 
theory models with a gravity dual. In section 5, we argue that using the density of states 
for a typical confining gauge theory with gravity dual in our models generically results in 
behavior qualitatively similar to the expected field theory description of an evaporating black 
hole, even though our model incorporates very few features of a realistic field theory. We 
consider a specific model and employ PDE techniques and numerics to verify our qualitative 
expectations for the evolution of energy distributions within the model. 

2 Basic setup 

At the most basic level, a field theory is a collection of degrees of freedom arranged in a 
translationally-invariant way in some number of dimensions. These degrees of freedom can 
be excited to different configurations depending on how much energy is available to them. 
As a simple model of this, we consider a lattice of sites, with some unspecified degrees of 
freedom on each site that can exist in p{n) possible configurations given energy n (which 
we take to be integer- valued for now). The configuration of the model at a given time is 
specified by the collection {{rir, K)} where n^ is number of units of energy residing on site 
r, and kr < p{nr) is a positive integer specifying the state of the degrees of freedom on 
site r. To incorporate dynamics into the model, we need to specify some transition rates 
{{ur, kr)t} -^ {{nr, kr)t+dt}- Wc constrain these by locality, conservation of energy and the 
principle of detailed balance, but otherwise simply postulate that all transitions consistent 
with these basic constraints are equally likely. Our construction gives rise to a particular 
set of stochastic equations for the evolution of a probability distribution of p({^r}) over 
the possible energy distributions in our model. We find that the most general evolution 



equations satisfying our constraints depend on only two detailed features of the underlying 
physics: the density of states p{E) for each site, and another function C{E) related to the 
energy-dependence of transition rates. 

2.1 Master equations for evolution 

To define a precise stochastic model using our basic setup, we assign a probability p{{nr}) 
to a particular energy distribution {ur}. We do not keep track of the detailed state, labelled 
by kr, of each site. The time evolution of such an ensemble is governed by the transition 
rates W^{{nr}^{n;,})5 giving the probability per unit time that the system in configuration {n^} 
will make a transition to configuration {n[.}. The time evolution is governed by a master 
equations 

dtP{{nr}) = Y^ [p({<})W^{<}^{n,} -p({nr})H^{„,}->{<}] • 
{<} 

For our system, we assume energy conservation and nearest neighbor interactions, so that 
the rate vanishes unless the total energy is conserved and the final state differs from the 
initial state only at a neighboring pair of sites. We assume that locality in the interactions 
further implies that the rate for a transition {. . . ,na,nb, ■ ■ ■} -^ {. . . , n'^, n'^, . . . } involving 
neighboring sites a and b does not depend on the state of the degrees of freedom on the 
other sites. Finally, we assume a discrete translation/rotational invariance such that this rate 
W(na,n,,)^{n'^,n') dcpcuds ouly ou na,nb,n'^, and ra'^, and not on the specific sites a and b 
involved. With these restrictions, the master equation takes the form 

dtPiirir}) = YY1 [^({^^1^'^ + k,nb- k})W^n^+k,n,-k)^{n^,n,) - pi{nr})W(^n,,n,)^(^n^+k,n,-k)] , 
{a,b) k^O 

(1) 

where {a,b) indicates nearest-neighbor sites and p{{nr\na + k,nb — k}) is the probability 
for the distribution that is obtained from {n.,,} by adding k units of energy to site a and 
subtracting k units of energy from site b. From this, we can give a general result for the 
evolution of the energy expectation value at a particular site 

Uc = Ypi{nr})nc . 

{rir} 

We find 

dtflc = ( X] X] ^ W(n,,na)^[n,+k,na-k) ) , (2) 

\{c,a) ky^O I 

where the first sum is over sites a that are nearest neighbors to c. 



^^For a general review of stochastic processes and master equations, see for example ifTSl . 



Equilibrium and detailed balance 

The master equations above act independently in each sector with a given total energy. In 
equilibrium, a basic assumption of thermodynamics is that all states with a given energy 
are equally likely. If p{n) is the number of states on a given site with energy n, then the 
number of states with the distribution {n^} is YlrPi^r), so we should have for an equilibrium 
configuration 

p({nj) =po JJpK) (3) 

r 

where po is a constant that can depend on the total energy. For such a configuration, the 
right-hand side of (H) must vanish, but we should generally expect a stronger condition, the 
Principle of Detailed Balance. This requires that in a given time, the number of transitions 
in the ensemble from distributions {n^} to {n[.} should equal the number of transitions from 
{n'^.} to {ur}. This is equivalent to saying that for the equilibrium configuration (|3]), each 
term in the square brackets in dD should vanish separately, i.e. 

Using (|3]), this gives a condition on the transition rates, that 

p{na + k)p{nb - k)W(na+k,n,-k)^(na,n,) = pina)pinb)W(na,ni)^{na+k,n,-k) ■ (4) 

Transition rate democracy 

Since we are interested in effects that do not depend strongly on the underlying microscopic 
dynamics, we will simply assume that the transition rate from a given initial energy distri- 
bution (ria, rife) of a pair of sites to various possible final distributions (ua + k.Ub — k) is 
proportional to the number of available final states with that distribution. This gives: 

W^K,ni,HK+fc,n6-fc) = C{na, nb)p{na + k)p{nb - k)f{k) , (5) 

where we allow a function f(k) (taken to be even) that can be used to introduce an additional 
dependence on the amount of energy transferred. As an example, we might want to allow 
only transitions up to a certain maximum energy transfer, using a function 

f{k) = 9{E^,^-\k\) 

Combining this with the constraint dH) above, we find that 

C{na, Ub) = C{na + k,nb- k) f{k) ^ 

but using this relation repeatedly givesG 

C{na, Ub) = C{{na + nb)/2) , 

where we have included the factor of 1/2 for later convenience. Thus, we have 

W^ria,ni)^{na+k,ni-k) = C((?^a + nb)/2)p{na + k) p^Ub - k)f{k) . (6) 



"^Here, we are assuming that any configurations with the same energy can be connected by a series of 
allowed transitions i.e. that there are no superselection sectors. 



Choices for W 

We now comment on a few specific choices for transition rates of the form ^. The simplest 
choice is 

W(na,nt)Mna+k,ni^k) = C p^Ua + k) p^Uh - k) . 

This model arises if we simply assume that all microscopic transition rates between the 
individual quantum states states are equal. 

We can also restrict the number of units of energy transferred in a step to obtain 

W^na,nb)^{na+k,nt-k) = CO{E^^y, - \k\)p{na + k)p{nb - k) . (7) 

For another physically motivated solution, consider the situation where nearest neighbor 
pairs interact with some rate q (independent of the pair and the energy on the sites) and 
where in each interaction, the two members of the pair instantly thermalize. In this case, the 
state after the interaction will be (n'^, n'^) with probability proportional to p{n'^p{n'^) with 
the probabilities summing to one. This gives a transition rate 

..r _ ^ p{na + k)p{nb-k) 

which is also of the form (|6]). With a constraint on the number of units of energy transferred, 
this gives 

Physically, the denominators in these expressions arise because we have introduced two time- 
scales into the problem; the final rate is the product of a basic reaction rate with a probability 
for various outcomes of the reaction (because the interaction that determines the final state 
happens very fast such that the probability for changing to a different state assuming the 
pairs interact is not much smaller than one). 

We can obtain a simpler model that should have similar qualitative behavior to ^ by 
replacing ^^ p{na + l)p{nb — I) with p^{{na + rib)/ 2). This gives 



W^K,n6H(na+fc,n6-fc) = C9{E^i,^ - \k\)^-^^ , ^ WON " ^^^ 



, p{na + k)p{nb - k ) 
p^iiria + nb)/2) 



2.2 Static distributions for Hagedorn density 

We would now like to analyze the physics of the class of master equations derived in the 
previous section. Starting from the general form (|6) of the transition rate, the equation (|2]l 
for the evolution of the average energy at a site becomes: 

9in,= /^J]A;C(K + r2j/2)p(n, + A;)pK-A;)/(A;)\ . (10) 

\(c,a) ky^O I 



In the next section, we will study this equation in the thermodynamic limit where fluctuations 
are small and derive a differential equation for the evolution of the energy distribution in a 
continuum limit where variations are small on the scale of the lattice spacing. However, 
we note that even at this stage, a special role is played by the Hagedom density of states 
p(n) = Be"^^. With this choice, we find 

\{c,a) k^O 
\(c,a) fc^O 

= 0, 

where in the last line we have used the fact that / is an even function of k. Thus, for an 
exact Hagedom density of states, the average energy at each site is constant. This is simple 
to understand: if we have two neighboring sites, the number of states with energy n^ on 
one site and n^ on the other depends only on Ua + Uh for a Hagedom density, so by our 
assumptions above, it is equally likely for a unit of energy to move from a to 6 as it is for the 
energy to move from 6 to a. 

For the result we have just derived to hold exactly, the Hagedom density must hold even 
for negative energies; more realistically, we should restrict to positive energy i.e. we should 
take p(n) = for ri < 0. In this case, starting with some distribution for which all sites have 
nonzero energy, the system will evolve with no change in the average energies until some 
members of the ensemble have zero energy for one or more sites. Since we cannot transfer 
any more energy from the sites with zero energy, these sites will gain energy on average, so 
there will be a slow diffusion of energy to sites for which the ensemble includes states with 
zero energy on these sites. 

3 Macroscopic equations in the continuum limit 

Retuming to an unspecified density of states, we will now analyze the equation ^ for the 
evolution of the average energy in the model © with transition rates of the form (|6]l. We 
have 

5in,= /^^A;C(K + nJ/2)pK + fc)pK-A:)/(A;)\ . (11) 

\(c,a) ky^O I 

For the large A^ theories that we consider (or more general theories where the individual sites 
represent enough degrees of freedom so that a thermodynamic limit can be assumed), we can 
replace the expectation value of the function of n^ here with the function of the expectation 
value to get an evolution equation for distribution of average energies, 

9t^c = J2J2^ ^((^-^ + ^a)/2)pK + k)p{na - k)fik) , (12) 

(c,a) k^O 

where now all ns represent ensemble averages. 



Let us now take a continuum limit of this equation. Specializing to one spatial dimension, 
we can take c = x and a e {x + e, x — e} to obtain 

dtn{x) = y^ kf{k) {C{{n{x) + n{x + e))/2){p{n{x) + k)p{n{x + e) - k) 

fc^O 

+ C((n(x) + n{x - e))/2){p{n{x) + k)p{n{x - e) - k)} , 

We also rescale energies k ^ 5k so that k becomes continuous in the limit 5 — )■ 0, and 
assume that f{k) falls off fast enough that only infinitesimal values of k contribute in the 
limit. In this case, taking the limits 5, e — )■ while rescaling C to leave a finite result gives 

dtE{x) = d, {-p\Eix))C{E{x))d/^''^P^ ^ 



dE j 

where we have denoted the continuum energy distribution by E(x). Generalizing to higher 
dimensions, we find 

dtE{x) = V ■ \^-p'{E{x))C{E{x))V^^^ . 

It is interesting that the right side vanishes identically if and only if we have a Hagedom 

density of states p{E) oc e^^. 

The expression \n{p{E)) is (up to a constant) the microcanonical definition of entropy, 

so we have 

d , , ,^,, dS 1 
J^ln(pm) = - = - = ^. 

In the end we have simply 

dtE{x) = V ■ {-p^iE{x))C{E{x))VP} . (13) 

For the simple model (JT]), we obtain 

dtE{x) = CV ■ {-p^{E{x))VP{E{x))} . (14) 

while for the model (|9]), we get simply 

dtE{x) = -CV^/3{E{x)) . (15) 

We note that by defining /3 with the property that 

^'{E) = p\E{x))C{E{x)W{E) 

the general model (fT3]) reduces to 

dtE{x) = -V^P{E{x)) , (16) 

the same form as the simplest model (fTSl ). Since p and C are positive functions, $'{E) has 
the same sign as (3'{E) for any E. We will see that the qualitative dynamics in the model is 
largely determined by the sign of /3'{E) as a function of E, so studying the simplest model 
with /3 = /3 should give us a good idea of the general behavior. 

8 



Perturbations around a uniform density 

As a start towards understanding the qualitative behavior of our macroscopic equation, con- 
sider small perturbations about a uniform energy distribution (which is a time-independent 
configuration for any density of states). Expanding (fT3T ) around a uniform density Eq, we 
find that the equation governing perturbations is 

where e = E — Eq. Thus, we get exactly the heat equation or the inverse heat equation de- 
pending on whether d^ ln(p)/d£'^ = rf/J/rfi? is negative or positive, with no time dependence 
for the Hagedom density. In summary, we have the following possibilities: 

• d(3/dE < 0: diffusive behavior, inhomogeneities decrease over time 

• d(3/dE > 0: unstable clustering of energy, inhomogeneities increase overtime 

• dfi/dE = 0: static energy distribution 

In realistic field theories, the density of states may be a function of energy that exhibits more 
than one of these behaviors. In this case, we can have different qualitative behaviors for the 
evolution depending on the local energy density. 

4 Density of states in large A^ field theories 

We have seen that Hagedom densities of states play a special role in our simple models, 
leading to static energy distributions, or equivalently a trapping of energy in certain subsets 
of degrees of freedom. This behavior is reminiscent of the expected behavior of black holes 
in gravitational theories. There, large amounts of energy sent in from an asymptotic region 
can form a black hole; the energy is trapped in a localized region (presumably described 
by a subset of the degrees of freedom in the fundamental description) for a time scale that 
goes to infinity in the classical limit. In this section, we recall that Hagedom densities of 
states appear naturally in large N gauge theories, some of which provide non-perturbative 
descriptions of quantum gravitational theories via the AdS/CFT correspondence. Thus, the 
Hagedom trapping of energy in these theories may be part of the explanation for the existence 
of black hole states in the dual gravitational theories. 

Consider first a simple model with two matrices X and Y of degrees of freedom, with 
each matrix element describing a harmonic oscillator [72]. We assume no interactions, but 
demand that the states are invariant under gauge transformations that act on the matrices as 
X — !■ UXU^^ and Y — t- UYU^^. If we define matrix creation operators A]^ and Ay in the 
usual way, then the gauge-invariant states are constructed with traces of products of these 
operators acting on a Fock space vacuum state. If each oscillator carries a unit of energy, 
then the number of single-trace states with energy E is on the order of 2^ = e^^^^'^\ since 
we must have N operators in the trace and each can be A'-^ or Ay. Some of these states are 
equivalent by cyclicity of the trace, but the overcounting is by at most a factor of E, so the 



leading asymptotic behavior of the density of states is Hagedom. Including multiple trace 
states does not affect this conclusion. 

As discussed in ifTTlfTll . this Hagedom behavior is obtained in the weak coupling limit 
of any large A^ gauge theory with a discrete set of matrix degrees of freedom (e.g. the 
momentum modes in a finite volume field theory). However, the analysis breaks down for 
energies of order N'^. For traces with of order N"^ matrices, algebraic identities allow us to 
rewrite certain traces in terms of other traces (e.g. Tr(M^) = |Tr(M)Tr(M^) — iTr(M)^ 
for a 2 X 2 matrix) and the result is that the density of states for energies of order A^^ crosses 
over to e^" where a < 1. 

Away from weak coupling, the expected behavior of the density of states in large A^ 
theories with a discrete spectrum of states (for example, large N CFTs compactified on a 
sphere or a circle) was discussed in detail in [[T2l (see sections 3.1, 6.5, and 7.2) and in [[TSl . 
appendix A. 2. Here, the behavior may be inferred from the gravitational dual description. 
Typically we have \r\.{p{E)) ~ E^ with a < 1 for energies of order 1 (corresponding to ex- 
citing perturbative gravity / field theory modes in the bulk), then a = 1 (Hagedom behavior) 
for energies of order A 3 where A is the t'Hooft coupling (corresponding to exciting string 
excitations in the bulk). At energies of order N"^, we first encounter a regime with a > I 
(corresponding to small black holes in the bulk) and finally, above a certain energy of order 
iV^, this behavior crosses over to a = d/{d + 1) (corresponding to large AdS black holes). 

Thus, for large N gauge theories with a mass gap, at both strong and weak coupling, we 
have a Hagedom regime in the density of states at intermediate energies, with slower growth, 
p{E) ~ exp(i?"^^), for energy densities of order N"^ and for energy densities smaller than 
order A* in strongly coupled theories with gravity duals. In the next section, we will explore 
the dynamics of energy distributions for such densities of states in our simple model. We 
focus on spherically symmetric energy distributions which initially have a large density of 
energy in a small-volume region. We will see that the existence of an intermediate-energy 
Hagedom regime in the density of states dramatically increases the time for the energy to 
leak away to infinity. 

This is in accord with our expectations for the gravity dual of such a theory, where a 
very localized collection of energy should correspond to a bulk configuration that will form 
a black hole. For gauge theories with a mass gap, the bulk dual has an infrared wall to which 
the black hole will fall and then very slowly evaporate. 

Note that for conformal holographic field theories on R'^'^, the bulk geometry dual to the 
vacuum state is the Poincare patch of AdS space. In this geometry, massive objects fall into 
a horizon in finite proper time that is typically much shorter than the black hole evaporation 
time. Since only this shorter time scale is visible in the field theory evolution, it should be 
difficult to deduce the long lifetime of gravitational bound states from the field theory in 
this conformal case. Here, the dynamics of an energy distribution would be expected to be 
dominated by a simple spreading behavior that corresponds to the bulk object falling towards 
the horizon. These expectations are consistent with the fact that the entropy density for a 
conformal field theory is related to the energy density (according to conformal invariance) 
by S* ~ i^;'^/(''+i). This suggests a density of states p{E) ~ exp(_E''^/*^'^+^') that corresponds 
in our model to a simple diffusive (spreading) behavior. 
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5 Evolution of energy distributions 

In this section, we study the evolution of simple energy distributions in our model, choosing 
the density of states to be similar to those in large N gauge theories with gravity duals that 
support long-lived evaporating black hole states. 

As we have shown above, our most general model (fTST ) can be rewritten to take the form 

dtE{x) = -V^J3{E{x)) , (17) 

where /3 is a function of energy for which $'{E) has the same sign as the derivative 13 {E) = 
dS/dE. 

We begin by comparing (fTTT ) with the energy conservation equation 

dtE = -V-Je. 

We see that the physical content of (fTTI) may be summarized by saying that the energy current 
is 

Je = V0iEix))) = P'{E{x))VE{x) . (18) 

For typical field theories, 13' {E) < 0, so energy flows from regions of higher energy to 
regions of lower energy. This is standard diffusive behavior. 

In theories with a regime in the density of states where (3'{E) > 0, (fTSi ) shows that the 
energy flows from regions of lower energy to regions of higher energy. Thus, with a very 
rapidly growing density of states, we expect a clustering behavior in which local regions of 
larger energy draw even more energy from the surrounding regionsO As we discussed in the 
previous section, this behavior in the density of states is found in holographic theories whose 
dual gravity description includes small black hole states. Our observations suggest that the 
(3'{E) > behavior may in fact underlie the dynamics that allows the formation of these 
small black holes. 

For theories with a Hagedom regime in the density of states, such as the large N theories 
discussed in the previous section, (fTSl ) shows that the energy current is zero in regions where 
the local energy density is in this Hagedom regime. In this case, dynamics is frozen in the 
regions of space where the energy density corresponds to Hagedom behavior. We will see 
below that the boundaries of such regions can move, so energy is not trapped forever as long 
as there is a lower-energy regime in the density of states with slower than Hagedom growth. 

Our observations here are consistent with those at the end of section 3 based on the 
analysis of small perturbations around a uniform energy distribution. 

5.1 Spherically symmetric initial energy distributions 

In this section, we consider the evolution of smooth spherically symmetric energy distribu- 
tions which are monotonically decreasing in the radial direction. 



^This can be understood very simply. As an example, consider a situation with energy E on each of two 
sites. If the number of possible configurations for 2E units of energy on a single site is larger than the square 
of the number of configurations for E unit of energy on a site, then by our democratic evolution rules, the 
system is more likely to end up in the configuration with 2E units of energy on a single site. The condition 
p{2E) > p{E)'^ is equivalent to saying that the growth of log p is faster than linear 
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Figure 1 : Energy dependence of (3 for density of states with a Hagedom regime for E > i^n- 




Figure 2: Evolution of an energy distribution with a Hagedorn regime in the density of states 
for E > Ell. The interface between Hagedom and non-Hagedom energy density moves 
inward with time, allowing the energy to escape eventually. 

For a density of states with (3'{E) < at all energy scales, the energy current points 
radially outward everywhere, and the energy diffuses away to infinity. If (3'{E) is bounded 
above by some negative constant —C, then we have 

\Je\ = \^'{E)\\VE\ >C\VE\ 

so the magnitude of the energy current is at least as large as for the ordinary diffusion equa- 
tion with diffusion constant C. Thus, an initially localized collection of energy will spread 
out at least as quickly as a system governed by some diffusion equation. In particular, if a 
certain amount of energy Eq is contained in a region of size Rq, the amount of time before 
the energy in this region is eEo can be bounded above by a time that is independent of Eq. 
This follows because the diffusion equation is linear, so if / is a solution of the diffusion 
equation then so is Eq/. Thus the timescale in /, which is set by the spatial profile (and in 
particular the scale Rq) and the diffusion constant, is independent of the overall energy E^. 

Effects of a high-energy Hagedorn regime 

Now let us consider the effects of a Hagedom regime in the density of states at energies 
E > Ell, with /3'{E) < for smaller energies, as depicted in figure ([I}- In this case, consider 
an initial energy distribution for which the central energy is larger than En, as shown in 
figure |2l The energy current is strictly zero in the regions with E > En, while energy flows 
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outward in the surrounding region where E > Eu- One might expect that the energy in the 
Hagedom region is trapped there forever; however, integrating (fTTI) over a ball with radius R, 
we find that flux of energy going out of the ball has a finite limit as R approaches the radius 
i?H where E = E^i. Thus, the interface point must move inward with time, as indicated by 
the dashed arrow in figure |2] The resulting energy distribution necessarily exhibits a step- 
function behavior at R^, since the energy distribution for R < R^ must be exactly the same 
as the initial energy distribution. 

As an explicit example, we consider in Appendix A the time evolution of a one-dimensional 
energy distribution that is initially 



E{x,0) 



Eq \x\ < a 

Ixl > a ' 



comparing the evolution with constant /3'{E) = —C (which gives a standard heat equation) 
to the evolution with (3'{E) = —C below energy En but l3'{E) = (a Hagedom regime) 
above energy En- We are able to find an exact solution for the evolution in both cases. 
Comparing the energy in the region |x| > a as a function of time for the two cases, we find 
that when Eq ^ En, the behavior for early times is 



2Eq\ - no Hagedorn 
E{\x\ >a,t) = { V^ 

AEu\/- with Hagedorn 



Thus, the rate of energy flow away from the region where the energy is initially concentrated 
is parametrically smaller when the density of states exhibits Hagedom behavior for energies 
above a threshold. 

Effects of a low-energy Hagedorn regime 

Next, we consider the effects of a Hagedom regime for low energies E < En, with (3'{E) < 
for higher energies. In this case, for an initial energy distribution with central energy 
larger than En, the energy current is strictly zero in the regions with E < En, while energy 
flows outward in the central region where E > En- Since the outward flux of energy just 
inside the radius Rn where E = En is finite, the interface point must move outward with 
time, allowing the central concentration of energy to spread eventually to a uniform disc 
with energy density En, as shown in figure |3l As above, the energy distribution necessarily 
exhibits a step-function behavior at Rn, since the energy distribution for R > Rn must be 
exactly the same as the initial energy distribution. 

Effects of an intermediate-energy Hagedorn regime 

Finally, we consider the effects of an intermediate-energy Hagedom regime, with /3'{E) = 
for En < E < E-p and 13' {E) < for both high and low energies. This behavior of /3 is 
similar to what we expect for certain holographic confining gauge theories, as we discussed 
in section 4. In this situation, we have a combination of the behavior in the previous two 



13 




E» 



J = 



Figure 3: Evolution of an energy distribution with a Hagedorn regime in the density of states 
for E < Ell. The interface between Hagedorn and non-Hagedom energy density moves 
outward with time, resulting in a flattened distribution with energy density En. 




Figure 4: Evolution of an energy distribution with a Hagedorn regime in the density of states 
for Eii< E < Ep. 

cases, as shown in figure HI The energy spreads outward in the central region and in the outer 
region, requiring the inner boundary of the Hagedorn region to move outward while the outer 
boundary of the Hagedorn region moves inward. These eventually come together, leaving a 
discontinuity at which the energy jumps from En to Ep. We expect that the central distri- 
bution eventually spreads out to one with constant energy density Ep, after which it evolves 
as we discussed in the case of a high-energy Hagedorn regime. To check this intuition, we 
have performed numerical simulations of the evolution equation in our model for (3{E) with 
this behavior. The results, described in detail in section 5.2, confirm the qualitative picture 
in figure |4l 

Comparison with expectations from holographic confining gauge theory 

The behaviour shown in figured is qualitatively similar to what we expect for the evolution 
of a highly concentrated energy distribution in a holographic confining gauge theory. Such 
an initial energy distribution might arise from a high-energy collision of glueballs. In this 
case, we expect that the large initial concentration of energy spreads out relatively quickly 
to a metastable "plasma ball" IfTSl with uniform energy density. At this stage, the energy 
distribution would be similar to the third frame in figure S] The plasma ball then slowly 
evaporates, with the boundary of the deconfined plasma region moving inward until the 
plasma ball has completely evaporated. 

The formation of a plasma ball and its subsequent evaporation corresponds in the dual 
gravity picture to the formation and evaporation of a black hole sitting at the infrared wall of 
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Figure 5: An inverse temperature function with E on the horizontal axis and (3 on the vertical 
one. It shows smooth transitions into and out of the Hagedom phase and was plotted for 

Ef = 20. 

the dual spacetime. Thus, the metastable black hole in the gravity picture corresponds in our 
model to the uniform energy distribution in the central region in the third frame of figure |4] 
In our toy model, we have seen that the Hagedom behavior (3'{E) = for E > En leads to a 
parametric suppression of the rate of energy loss from such a configuration (at least at early 
times). Hence, it is plausible that the Hagedom regime in the density of states for the actual 
holographic confining gauge theory may be part of the explanation for the long lifetime of 
the black hole state in dual gravitational theory. 



5.2 Numerical results 

In this section, we present some numerical results for the evolution of spherically-symmetric 
energy distributions in a model for which the density of states has an intermediate-energy 
regime with nearly Hagedom behavior, in order to check the qualitative expectations dis- 
cussed in the previous section. 

It is possible to treat the Cauchy problem for the equation (fTTI) with finite difference 
schemes. We have used the Crank-Nicolson method, a second-order approximation stable 
for large timesteps, to demonstrate some of the qualitative features described in this paper 
In its most well known form, Crank-Nicolson is used to simulate the heat equation. Since 
it is an implicit Runge-Kutta method, it requires a system of linear algebraic equations to 
be solved at each time. Applying the Crank-Nicolson method to our equation results in a 
set of non-linear algebraic equations which we solve by Newton's method. Figure |5] shows 
the /3{E) function used in these simulations. Motivated by the behavior for holographic 
confining gauge theories, we choose /3 to decay with a power of — ^i^ at low energies and 
with a power of — ^ for high energies E > Ep with a nearly flat region in between. We do 
not include a regime with /3' > 0, since this complicates the numerics significantly. 
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Figure 6: Different initial conditions both evolve towards a state that has a flat line of energy 
at Ep before decaying at the edges. In the special case where (3{E) is a pure power law, the 
observation of all initial conditions approaching the same limit is known to hold precisely. 
Logarithmic scales are used on both axes in order to make the important features visible. 

Most of the initial conditions that were tested had the form 



E{x,0) = 20Ef 



1 + x^ 



chosen to give a smooth, localized energy distribution with central density significantly 
greater than the value Ep corresponding to the crossover to the high-energy regime in the 
density of states. 

We present the results of a few of our simulations in figures |6]and|7l 

The simulations we have performed exhibit the following general behavior: 

• At early times, the evolution depends in detail on the initial condition. However the 
distributions generally approach one with an approximately uniform energy density 
Ep in a central region, separated by a sharp boundary from a region of low energy. 

• The evolution from the initial energy distribution to an energy distribution with a 
plateau at E r^ Ep happens on a much shorter timescale than the subsequent decay of 
this flattened distribution. 

These observations are consistent with the qualitative picture described in the previous sec- 
tion and summarized in figure H] In particular, the last observation matches our expectations 
for holographic confining gauge theories (or their dual gravity theories) where the evolution 
of an initial energy distribution to a nearly-static plasma-ball (black hole) configuration is 
expected to occur much faster than the subsequent evaporation of the plasma ball (black 
hole). 

The discussion in this section has been mainly qualitative, backed up by numerical sim- 
ulations. In appendix B, we include a discussion of some rigorous analytic results pertaining 
to solutions of general equations of the form (fTTj ). These may be useful in future investiga- 
tions. 
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E(x.i) 




Figure 7: Evolution of an energy distribution for (3{E) similar to that in figure |5] The 
formation of a sharp boundary between high and low energy regions is a robust feature that 
appears for a wide range of initial conditions. 

6 Discussion 

In this note, we have illustrated by use of a simple stochastic toy model a connection between 
the density of states for subsystems of a system with a network of locally interacting degrees 
of freedom and the qualitative dynamics of energy distributions in the system. Depending 
on whether the density of states for sites in our model has a growth rate slower, faster, or 
equal to Hagedom growth {p{E) ~ e^^) we find diffusive, clustering, or static behavior for 
the energy distribution. 

It is an intriguing possibility that our observations may help us understand why certain 
large N gauge theories are able to describe the physics of black holes in their dual gravi- 
tational description. One of the unusual features that distinguishes large N gauge theories 
from more typical field theories is the existence of Hagedom and super-Hagedom regimes in 
the density of states. We have seen, within the context of our toy model, a direct connection 
between these behaviors and trapping / clustering of energy. Thus, it is plausible that in the 
full large N field theories with gravitational duals, the formation and long lifetime of black 
hole states may also be at least partly explained by the behavior of the density of states for 
subsystems of the field theory. 

While our toy model appears to give results for the dynamics of localized energy dis- 
tributions that are qualitatively similar to our expectations from certain field theories, we 
emphasize that this type of model is far too simplistic to give an accurate model of the field 
theory more generally. For example, while we have incorporated conservation of energy 
into our model, we have not included momentum as a conserved quantity. Thus, for ex- 
ample, the toy model cannot accurately model radiative transfer of energy. In the future, it 
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may be interesting to consider more detailed models that incorporate additional features of 
the field theory to see whether we can find behavior matching a larger class of gravitational 
phenomena. 

It is interesting to note that the macroscopic equations (generalizing (fTTT l) for a more 
general stochastic model including conserved momentum would take a form similar to hy- 
drodynamics equations. At least for near-equilibrium thermal states (and perhaps more gen- 
erally for large A^ theories) there is a general expectation that hydrodynamics should indeed 
provide an accurate description of the field theory dynamics. Furthermore, in recent work, 
specific hydrodynamical models associated with large A^ gauge theories have been shown 
to have solutions in one-to-one correspondence with solutions of Einstein's equations ex- 
panded around black brane solutions in gravity with a negative cosmological constant [[T6l . 
It would be interesting to see if one could reproduce these hydrodynamics equations as the 
macroscopic equations of a stochastic model for which the transition rates depend on both 
momentum and energy. In this case, one might be able to understand more directly how 
the hydrodynamic coefficients (viscosity, etc..) arise from microscopic properties of the 
underlying field theory. Additionally, such a construction would demonstrate that models 
only slightly more complicated that the one we considered here are capable of giving rise to 
detailed gravitational dynamics. 
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A Example of diffusion with a Hagedorn regime 

In this section, we consider a simple example to see the timescale for diffusion of a localized 
collection of energy in a case where the ordinary diffusion equation governs diffusion for 
E < Eu while the diffusion constant drops to zero (as in a Hagedorn regime for our model) 
for E > Ell. Thus, we would like to find the solution of 

dtE = dl^{E) 

where 

^^E) = 1 ^ E<Ei, 

with initial conditions 

^^^'°) = \i?o \x\<a ' 

where we assume that E^ is greater than E^. We would like to compare our result to the 
simple diffusion case with $(£') = E m order to see the effects of having a Hagedorn 
regime. 
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As discussed in section 5, the solution for t < T is expected to be discontinuous at 
some points ±RE_{t), with E = Eq for \x\ < Ruit) and E decreasing from E = E^ as 
|x| increases from Rn. The point of discontinuity -Rh(^) decreases with time from x = a, 
reaching x = at time T. 

Suppose at some time we have -Rh(^) = R- Since the energy density jumps from Eq 
down to Ell at a: = Ru, the rate of change of the energy in the region < x < i? at this 
instant will be —{Eq — £'h)-Rh(^)- This must equal the energy current just to the right of Ru, 
so we must have 

{Eq - EH)Mt) = d.,E{R^, t) . (19) 

Our strategy will be to find a solution e{x,t) to the ordinary diffusion equation such that 
the point Ru{t) at which e{Rii{t),t) = En satisfies the equation (fT9l ). In this case, the 

function 

e(|x|,t) |x| > Rn{t) 



E{x,t) — A rn II n /aN 

^ ' \ Eq \x\< RH{t) 

will be a solution of our original differential equation, since it is static for x < Ru, solves 
the ordinary diffusion equation for x > R^, and satisfies the boundary condition (fT9] ) by 
construction. 

We will show that the desired function e{x, t) is the solution to the ordinary diffusion 
equation with initial conditions 

/ n\ f X > a 

^ ' ^ [A -cx) < X < a ' 

where A is a constant determined by Eq and Eu- For these initial conditions, we have 

e(x,r) = , / axe « 

A / ^ f a- X 

1 + err 



2 V V 2Vt 

Now, choosing A so that E^ is between A/ 2 and A, we can define Rn{t) as above by 

En = e{Ru{t),t) , 



which gives 



Now, using 



we can check that 



R^{t) = a-2^/t eir^ (2-^ - 1 j 



OxC = ; e 

Vivrt 



d,e{Rn{t),t) A e-(-f"(2EHM-i))2 

: = , f- = constant 

Ru{t) V^eTr\2En/A-l) 
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Thus, we obtain the desired solution by choosing A such that 

^ g-(erf-i(2BHM-l))2 

In summary, the solution is 



(Eo - En) = ,^ „.,,-!,, ^ ,, ,^ • (20) 



Eo \x\<Ruit) 



with 



R^(t) = a-2Vi err^ (2^ - 1 j 



and A given by (1201) . More explicitly, A is determined in terms of E^ and Eq by 



where / is the solution to 

v^/e^'(l + erf(/)) 



/2/T ■ --€rT\\ -^H 



-E'o — Ell 



For comparison, the solution of the ordinary diffusion equation with the same initial condi- 
tions is 

To see the quantitative effects of having a Hagedom regime for E > En, we can compare 
the flow of energy outside the region [—a, a] for the solutions ( 12T]) and ( l22l) . In each case, we 
calculate the energy in the region |x| > a as a function of time. With the Hagedom regime, 
we find that while Rii{t) > (i.e. while the energy density at a: = is still greater than En), 

~ 4EhW^ (^o>^h) 



For the simple diffusion case, we find 



= 2EoJ\l-e-'r) + 2Eoa(l-eii(^ 



2EoW ^ (t < a^) 
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Thus, at early times, the energy flow away from the central region is smaller by a factor 
2Eii/Eo with a Hagedom regime. 

We note that the time before the central maximum in the energy distribution decreases to 
the Hagedom threshold E^ (equal to the time T when -Rh(^) = 0) is 



1 / a \ 7T f Eq 



In the case of simple diffusion, the time when E{x = 0) = Eu is 

^- = (2erf-' (£„/£„)) " I ("f ) ■ <'"' 

This differs from (|23T ) only by a factor of order 1, so over longer time scales, the flow rate 
of energy out of the central region becomes comparable for the two cases. The reason is 
that as energy leaves the central region, it does not escape quickly to infinity as it would in 
a relativistic field theory, but rather diffuses slowly away. This results in a buildup of energy 
outside the central region which results in a smaller d^E at the boundaries of the central 
region, and thus a smaller flow rate. Thus, at these later times, the flow rate away from the 
central region is limited by the speed at which energy can diffuse away. This is the same in 
both cases, so the effects of the Hagedom regime in ^{E) is masked somewhat. 



B Concentration Comparison 

The equation 

dtE = dl^E) , (25) 

equivalent to (fTTI) . is commonly referred to as the Filtration equation and appears in the 
study of diffusion in porous media [[T8l . In this appendix, we point out a few general results 
pertaining to solutions of this equation that may be useful in future investigations. 

For smooth axisymmetric energy distributions that monotonically decrease in the radial 
direction, it follows from (fTTT ) that that the maximum energy density at r = decreases 
monotonically with timeO The qualitative behaviour is that energy diffuses from large den- 
sities at small radii, where V^-E < 0, to small densities at large r, where V^-E > 0. Note 
that (each component of) WE stays non-positive at all r ^ 0; therefore r = is the only 
local maximum at all times 



''This is a consequence of the so-called maximum principle. 

^For WE to become negative at ?■ 7^ it has to first vanish. The point where \/E\r^o first vanishes is an 
inflection point, i.e. the components of VV^-E are negative. Then, taking the gradient of the (fT6] l shows that 
WE remains non-positive: 

dtWE = -P'W^E < 0. 
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Consider a ball B(R) of radius R centered at r = 0. We denote by Mnit) the total 
amount of mass inside the ball: 



Jb(r) 



IB{R) 

This mass decreases as a function of time as: 

dtMn = -/3' / VE{t) < 0, (26) 

Js(R) 

where we have integrated (fTTT l and S{R) is the boundary of B{R). We can also define the 
diffusion time t^{R, e) to be the shortest time it takes for MR{t) to drop to e times the total 
mass, i.e. Mr{U) = eMoo(O). 

In order to obtain quantitative bounds on the time scales for evolution of energy distribu- 
tions for general ^(E), we can use the following comparison theorem [[T7l[T8l : 

Concentration comaprison. Let ei{x, t),i = 1,2,?, be solutions to dtCi = V^$j(ei), respec- 
tively, for monotonically increasing ^i{ei) and ^i{0) = 0. If ei{x,0) = e2(x, 0) = 63(0;, 0) 
and^i{e) h $2(e) h ^3{e) for all ^ then, Mi{R,t) < M2{R,t) < M3{R,t) for all Rand 
t, where Mi{R) = J^^^^ d'^x e^. Therefore: tl{R, e) < tl{R, e) < tl{R, e). 

As an example, we consider the evolution of an energy distribution in our model equation 
(fTTl ) in the case where (3{E) has Hagedom behavior (3'{E) = for E > E^ and typical field 
theory behavior (3{E) = E°' with < a < 1 for i? < E^i. In regions where the energy 
is less then E-^, the evolution equation (flTl) reduces to the well-known super-fast diffusion 
equation: 

dtE = a{l - a)V ■ (E'^-^VE) . (27) 



The super-fast diffusion takes its name from its singular behavior at i? — )• 0. From (1181) . we 
see that the energy current is 

Je = -a{l - a)E''-^VE, 

which diverges for E —^ 0. This singular behavior means that energy can escape to infinity 
in finite time. However, in realistic models, the density of states will be modified from the 
p{E) = exp{AE°') behavior at the lowest energies so that this behavior is avoided. Here, 
we will avoid the problematic low-energy behavior by considering initial Cauchy data with 
a small background energy density, 

E{x, t = 0) = e{x, t = 0) + Eb, 

where e(x, t) is a non-negative function that vanishes at spatial infinity. Therefore, the prob- 
lem we address here is the relaxation of far-from-equilibrium initial distributions evolving 



'./(e) t gie) if and only if /^ f{e')de' > J^ g{e')de'. 
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on a uniform background. Note that introducing a background, no matter how small, dra- 
matically changes the solutions to the fast-diffusion equation since it cuts off the propagation 
speed of perturbations 

In summary, we wish to consider a solution to the equation 



dtE = dl^{E) , (28) 



with 



with initial conditions 



HE) = 1 




1 E<En 
' E>E^ 


E(X,0): 


\ Eb 
\ Eo 


x\ > a 
x\ < a ' 


> ^H > Eb. 







where we assume that Eq 

If E(x,t) is the desired solution, then we note that e{x,t) = E{x,t) — Eb will be a 
solution to the equation (l28l) with $(e) = ^(e + Eb) + C, where C is an arbitrary constant. 
In order to apply the concentration comparison theorem, we will choose C so that $(0) = 0. 
Thus, we consider the equation (l28l) with 



and 



, , aE^~'-a{e + EB)^~' e < E^ - Eb ,.„^ 

^' ^ aE^~^-aE^^^ E>E^-Eb ^ ' 

Ixl > a 



e(x, 0) > „ rn I I 
^ ' ^ Eq- Eb \x\< a 



We would like to know how long it takes before the energy density is everywhere smaller 
than the Hagedom density. To obtain bounds on this time scale, we use the concentration 
comparison theorem, choosing $2(e) as in (|29l) , and taking $i and $3 to be of the form 



UE) 



BiE e<e; 

E> E. 



M 



M 



for which we obtained the solution in the previous section. By choosing Bi and Ej^'^ opti- 
mally so that the conditions of the theorem are satisfieco, we can show that the time before 
the energy density is less than E^ everywhere is bounded by 



:{aEo)' <T< ^—^{aE^) 



a)E^' ' - - AaE^ 



^In d dimensions there are no solutions to (l27T i for a < d/2 due to the phenomenon of instantaneous 
extinction iflTl . However, in the presence of the background, the solutions exist at least for all < a. 

'"Specifically, the optimal bounds come by ensuring that /g " ^ ^1 ~ lo" ^ ^2 = Jq" ^ ^3 and then 
choosing Bi and ^3 as large and as small as possible so that L $1 < /„ $2 < /n *&3 for all E < Eq — Eq. 
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